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Abstract 

We present numerical evidence that solutions of the Bethe Ansatz equations 
for a Bloch particle in an incommensurate magnetic field (Azbel-Hofstadter 
or AH model), consist of complexes — "strings". String solutions are well- 
known from integrable field theories. They become asymptotically exact in 
the thermodynamic limit. The string solutions for the AH model are exact in 
the incommensurate limit, where the flux through the unit cell is an irrational 
number in units of the elementary flux quantum. 

We introduce the notion of the integral spectral flow and conjecture a 
hierarchical tree for the problem. The hierarchical tree describes the topology 
of the singular continuous spectrum of the problem. We show that the string 
content of a state is determined uniquely by the rate of the spectral flow 
(Hall conductance) along the tree. We identify the Hall conductances with 
the set of Takahashi-Suzuki numbers (the set of dimensions of the irreducible 
representations of Uq{sl2) with definite parity). 

In this paper we consider the approximation of noninteracting strings. It 
provides the gap distribution function, the mean scaling dimension for the 
bandwidths and gives a very good approximation for some wave functions 
which even captures their multifractal properties. However, it misses the 
multifractal character of the spectrum. 

I. INTRODUCTION 

The problem of two-dimensional electrons in a periodic potential and a uniform magnetic 
field (Azbel-Hofstadter problem []I],0]) has a rich history and numerous applications. It is 
equivalent (in the Landau gauge) to a one-dimensional quasiperiodic difference equation: 

ijjn+l + V'n-1 + 2A COsiky + 2TTnT])lpn = El/jn (l) 
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with two competing incommensurate periods 1 and l/rj. Here rj is the flux of the magnetic 
fleld through the unit cell in units of the flux quantum h/e. This equation is also known as 
the Harper equation, the almost (or discrete) Mathieu equation, etc. It has been applied to 
quasicrystals, localization/delocalization transition quantum Hall effect and even 

DNA chains |^ . For recent review of models of Hofstadter type see [|] and references therein. 

The spectrum of this equation is complex. If the flux is a rational number P/Q, 
there is a common period Q in the equation (|ip and the spectrum consists of Q bands. 
In the incommensurate limit, when rj is an irrational number^ (P — > cxd, Q —>■ cxd) the 
spectrum becomes an inflnite Cantor setQ |T]J^ with total bandwidth (Lebesgue measure 
of the spectrum) 4|A — 1| [P,p!0|,pn]] . At |A| = 1 the spectrum becomes a purely singular 
continuous — uncountable but measure zero set of points fl^ (for review see |ll3|,|5||). There is 
numerical evidence that in this case the spectrum and wave functions are multifractal . 



The scaling description of multifractal sets, generated by a quasiperiodic equation, an 
important and challenging problem, is far from being understood. The reasons are twofold: 
a lack of language to describe scaling properties of singular continuous sets and a lack of 
methods to attack incommensurate systems. 

Recently, it has been shown that the Harper equation, belongs to the class of integrable 
models of quantum fleld theory. Despite being just a one particle problem, it has been 



"solved" by the Bethe Ansatz |T5[. This creates a perspective for an analytical solution of 
the problem and eventually for the employment of methods of conformal fleld theory for 
flnding scaling properties of multifractal sets. 

Below we make a flrst step towards solving the Bethe Ansatz equations for the most in- 
teresting critical case |A| = 10. We show that in the incommensurate limit roots of the Bethe 
Ansatz equations are grouped in complexes called strings, so that each state is a composition 
of strings. Strings are well known for standard integrable models of quantum fleld theory, 
for instance XXZ Heisenberg chain [T^ . They become exact only for a macroscopic system, 
where the number of particles and the size of the system are sent to inflnity. The role of 
thermodynamic limit for the AH problem is obviously the incommensurability — P, Q — > cxo, 
P/Q —>■ rj = irrational. The common period Q plays the role of the size of the system. 
Indeed, we present numerical evidence that the "string hypothesis" remains valid — strings 
become exact in the incommensurate limit Q —>■ oo. 

A flrst important application of the string structure of solutions is a detailed hierarchical 



tree of the spectrum, i.e. an algorithm for generating this Cantor set spectrum [18|. The 



^Although some properties of the spectrum depend on the type of irrational number rj, in this 
paper we consider the so called typically Diophantine numbers. These numbers have a full Lebesgue 
measure and thus suffice for almost all physical applications. 

^closed, nowhere dense set which has no isolated points 



^The initial progress toward solving Bethe equations of |15] has been made in where the 
explicit analytical solution for zero energy level as well as some numerical results for midband 
levels are reported. 
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hierarchical tree gives the topology of the set (Sec. We show that the string decom- 
position of a state is tied to the holonomy of the wave function, i.e. the Hall conductance 
of the state, and therefore must be of an algebraic geometrical nature. The set of Hall 
conductances generates the spectral flow, which in its turn describes the hierarchical tree 
(Sec. 0). 

Since strings become exact in the incommensurate limit with an accuracy 0{Q ^), at 
large but finite Q, the "bare" noninteracting strings give a very good numerical approxima- 
tion for some wave functions. This approximation even captures the multifractal properties 
of wave functions (Sec. [V|). The noninteracting strings give rise to a distribution function 
of gaps of the spectrum (Sec. |V A] ). The string hypothesis also gives a mean field scaling 
exponent for the spectrum (bandwidth distribution), and provides a framework for calculat- 
ing anomalous dimensions as "finite size" corrections to the bare values obtained at Q — oo 
from string solutions. We will not compute them in this paper. 

We start by the construction of the hierarchical tree for the spectral flow of a general 
quasiperiodic system (Sec. |ID. Then in Sec. |T| we review the Bethe Ansatz for the AH 
problem and formulate the string solution. In Sec. we argue that string solutions are tied 
to the Hall conductances and generate the hierarchical tree described in Sec. |I[ The string 
hypothesis is used to compute the gap and bandwidth distributions in (Sec. |V A| and Sec. 
|V B|) . We conclude with a list of results and "loose ends". In the Appendix (Sec. |^ we 
present a technically involved derivation of Bethe Ansatz equations. 



II. HIERARCHICAL-TREE: SPECTRAL FLOW, HALL CONDUCTANCE AND 

INTERMEDIATE FRACTIONS. 



The hierarchical tree is an adequate language to describe properties of quasiperiodic 
systems and strange sets these systems generate. In this section we describe the hierarchical 
tree for the spectrum of the AH problem, although a major part of the construction seems 
to be valid for more general quasiperiodic linear systems. Later in Sec. ^ we show that this 
tree determines the structure of the Bethe roots as well. 

The hierarchical tree may be defined as follows. Let us consider a sequence of rational 
approximants r]^^^ = P^^^Q^^^ with increasing Q^^^ to an irrational flux t] so that: \r]^^^ —v\< 
c{Q^^^)~'^, where c is j-independent constant^. We refer to the spectrum of the Harper 
equation characterized by the flux t]^^^ = P^^^Q^^^ as the generation of the hierarchy. Let 
us connect the fc-th band of the generation t]^^^ (the daughter generation) to a certain band 
k' of some previous (parent) generation r]^'^'\ thus creating an infinite graph with no loops. 
We call this graph a hierarchical tree if energies Ej{J) belonging to any branch J of the 
tree form a sequence converging to the point E{J') of the spectrum in such a way that the 
sequence {Q^^^)^-'J\Ej{J) - E{J) \ is bounded but does not converge to zero. The numbers 
e J are anomalous exponents^. The example of hierarchical tree for AH problem is shown in 



This sequence always exists and can be constructed e.g. from the so called Farey Series [21| 
^To avoid confusion let us stress that the very existence of the hierarchical tree is a (scaling) 
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Fig. |l]. To construct the hierarchical tree it is necessary to find the sequence of generations 
and a rule to determine the parent band out of a given band of a given generation. 

Below we conjecture that the hierarchical tree is the spectral hierarchy — an integral 
version of the spectral fiow. Let us recall the notion of spectral fiow and its rate — the Hall 
conductance. Consider the spectrum of the problem for a given fiux r/(^) = P/Q. It 
consists of Q bands separated by gaps. Let us choose some gap k and change the fiux by an 
infinitesimal drj. The spectrum remains unchanged except for some new levels that appear 
within a gap. The number 5Nk of levels crossing the energy lying inside the gap, is a 
spectral fiow. Its rate 

is known to be the Hall conductance of the gap [^. This number is an integer and depends 
on the gap. The difference between the Hall conductances of neighboring gaps, i.e. the 
spectral flow into the fc-th band a{k) = ak — Ck-i is the Hall conductance of the band. 

Let us now consider two close rational numbers i]^^^ = P/Q and ri^^~^^ = P'/Q' with 
Q' < Q. The number of states in a band is 1/Q and 1/Q' respectively in units of A^, where 
is the total number of lattice sites. If there is a band k of the problem with the flux P/Q, 
and we can flnd a flux P'/Q' with Q > Q' , P > P' such that its Hall conductance is a ratio 
between the difference of the number of states and the fluxes 

then we say that the band k of the generation P/Q has a "parent" band in the generation 
P'/Q'. This formula may be viewed as an integrated Streda formula @. It determines the 
flux P'/ Q' and by virtue of an iterative procedure, generates a sequence of rational approx- 
imants (generations) t]^^^ to the flux t]. The sequence t]^^^ differs from the approximants 
obtained by truncating the continued fraction expansion of 

V = , ^ 1 = [ni, 722,^3,...]. (4) 

^1 H ; — I — 

Instead it consists of so called intermediate fractions approximants (see below) ||23|| . 

We complete the integrated Streda formula by the adiabatic principle: let us enumerate 
all states from the bottom of the spectrum and characterize them by an integrated density 
of states (IDS) z/ = iVgtates/Q- Here Ai'states is a total number of states with energy less than 
the energy of the given state. All states in the k-th band have the IDS {k — l)/Q < v < k/Q. 



hypothesis. Moreover, the tree constructed in this paper is a conjecture. There is no deeper reason 
for the hypothesis and the conjecture other than physically plausible arguments and numerical 
evidence. 
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The adiabatic principle assumes that a middle stateQ in the band k with IDS {k — 1/2) /Q 
has a parent in the parent band k' such that {k' — 1)/Q' < {k — 1/2) /Q < k' /Q' . This gives 



k' 



Q^^ 2' 



+ 1 



(5) 



and together with integrated Streda formula determines the hierarchical tree. [x\ denotes 
the integer part (floor function) of x. 

To construct the tree, we need to know the Hall conductance of a band <j{k) and the 
Hall conductance of a gap au- The Hall conductivity of the fc-th gap at generation P/Q is 
known to be a solution of the Diophantine equation 1 2^,25] 



PcTfc = k (mod Q) 



(6) 



restricted to the range —Q/2 < < Q/2. 

To write the solution of this equation, let us represent the flux ri^^^ = P/Q of 
the given generation as a continued fraction (^ and denote the convergents Pi/Qi = 



[ni, n2, 
is 



Pr/Qr = P/Q 



V 



j) 



Then the Hall conductance of a gap 




(7) 



where {x} is the fractional part of x. 

In its turn the Hall conductivity carried by the k-th band cr{k) = — (Tk-i- It is 
easy to see that the Hall conductivity of the band takes only two values {—ly^^Qr-i and 
{—lyiQ — Qr-i)- It determines the parent generation by virtue of Q' = Q — |cr(A;)|: 



V 



P' 

0-1) = f_ 

Q' 



Pr Pr -- 
Qr Qr- 
Pr-l 



if a{k) 
if a{k) 



:-i)'-'Qr-i 

[-lyXQr-Qr- 



(8) 



This formula looks better in terms of the continued fraction 
Jj-i) _ j [ni,n2, . . . ,nr - 1] 



\ni,n2, 



, n. 



T-l 



if a{k) 
if a{k) 



ly-^Qr^i 

iy{Qr — Qr-l] 



(9) 



Thus the parent generation is obtained by either subtracting 1 from the last quotient of 
the continued fraction or by truncating the fraction. The sequence of generations produced 
by these iterations is known as the intermediate fractions . These numbers are deflned as 
numbers of the type n~^V^n~^ with u = 0,l,2,...,nj. The intermediate fractions list all 



the best rational approximants of the number rj p6|. This is sufficient to construct the 
hierarchical tree. 



^Middle state is the state for which the number of states in the given band with lower energy is 
the same as number of ones with higher energy. 
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Applying the same procedure to the parent generation we obtain another set of values 
of Hall conductances of bands. The entire set of Hall conductances generated by the flux 
P/ Q are numbers less than Q of the form 

{Qi-2 + mQi.i\m = l,...,ni, i = 1, 2, . . . , r}. (10) 

This set consists of numbers (denominators of intermediate fractions) known in integrable 
models related to Uq{SL2) as Takahashi- Suzuki numbers They are allowed lengths of 
strings in solutions, say for the XXZ model. These numbers are also the set of possible 
dimensions of irreducible highest weight representations of Uq{SL2) with definite parity 
27| . We might have expected these numbers to play a role, considering that the intimate 



connection between Uq{Sl2) and the AH-problem has already been outlined in previous 
pubhcations . 



Each path of the tree may also be characterized by a sequence of fractional numbers 
of bands: = k/Q^^\ lying on the path and converging to a given irrational fraction 
V. According to eq.(§) the parent fraction is determined by the daughter one as z/^-^"^^ = 
Q(h) {\Q^^~^\^^^^ - ^)] + l) with Q(^-i) = g(^) - \a{k)\. The sequence z/(^) converges 
to the irrational u faster than (Q*--'-*)^^, l^^'--'^ — u\ < c((5*--'^)~^. In terms of the fractions one 
may reformulate the scaling hypothesis as \Ej — E\ < c\u^^'^ — z/|"^"^\ defining the scaling 
exponent a{J). 

To illustrate the construction of the hierarchical tree described above, let us consider 
some particular flux, say, r] = P/Q = 4/15 = [3,1,3]. The spectrum of AH problem for 
this flux consists of 15 bands and 14 gaps. The ancestral generations of the tree are given 
by the sequence is ^ , |, |, |, f,^, Notice, that the truncated continuum fractions 
are different. They are y , | , \, — just a subset of generations. The Hall conductances of 
bands and gaps (ordered according to their energy) for each generation are: 



^ = — 

m = 

1 

V4 = - 

V3 = 
V2 = - 




-- 4,-11,4,4,4,-11,4,4,4,-11,4,4,4,-11,4 

- 4, -7, -3, 1, 5, -6, -2, 2, 6, -5, -1, 3, 7, 4 

- 4,-7,4,4,-7,4,-7,4,4,-7,4, 
-- 4,-3,1,5,-2,2,-5,-1,3,-4 

-3,4,-3,4,-3,4,-3 
-3,1,-2,2,-1,3 

1.1, -3,1 

1.2, -1 

1,-2,1 
1,-1 

1,-1 
1 



One should notice the asymmetry of a{k) for rj^. This asymmetry is related to the degeneracy 
present for the fluxes with even denominators (two bands touch each other at energy E = 0). 
The two possible choices ct2 = ±2 for the Hall conductance in the "gap" at = (actually 
there is no gap) restore the symmetry of the tree with respect to E —>■ —E. This tree, 
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depicted in Fig.|l| (half of the tree is shown) consists of 15 non-crossing branches Ji, i = 
1, 2, ... 15. They connect bands, characterized by the flux and number of the band, with 
Hall conductances (T{k) as listed below. The branches are listed in order of increasing energy. 
Only the lower half of the spectrum is presented 



i 1 1- i 1 1- i 1 1- - 

2? -^5 -^5 35 -^5 -^5 45 -^5 -^5 7 

J2 = (i,l,l;|,l,l;i,l,l;4,2,ll 



7 = ( L 1 I'll I'll — 1 A- — A A\ 



•-^3 (2'-'-'-'-i3'-'-'-'-i4'-'-'-'-ill'^' "^'IS'*^'^)' 

•^4 (2 3' 75 ^5 ^' 11 ' ^' 15 ' ^' n 1 "1 

? — I'i 9 — 9- i 9 1-2 Q _Q- A A A- A. ^ A\ ^ J 

•^5 V3' ' '^i4)'^;-'-)7)'-') ^iii;^5^ii5?'-'5^j) 

J& = (37 2, — 2; ^, 2, 1; 6, — 11), 

Ji — (3) 2, — 2; Y, 2, 1; 5, — 7; jg, 7, 4), 

/j, = fi 9 —2- - 4 4- — fi 4- — 8 41 

Here in brackets we showed the flux, the number of the band and the Hall conductance of 
this band. 



III. BETHE ANSATZ EQUATIONS AND THE STRING HYPOTHESIS 

The Bethe Ansatz solution is available for any point of the Brillouin zone < A;^. < 
27r/Q, Q < ky < 2ti {rj = P/Q) [Q, but it looks especially simple in the rational points of 
the Brillouin zone (the deflnition of rational points and corresponding values of kx and ky 
are given in Sec. [TX] ). They correspond to the centers and edges of bands. The study of 
these points is sufficient for our purposes^ Below we sketch the results of the Bethe- Ansatz 
solution. The details can be found in the Appendix and also in Refs. [T^,^^. The Harper 
eq. (|l|) at rational points can be mapped onto a functional equation 



fXKzE^(z) = iq{z + iTKq '^)(z — inq 2)vl>(g2;) 

—iq~^{z — iTK,q^){z + iKq^)'^{q~^z) (12) 



where q = e*'^^ and t, k, jj = ± are labels of the rational points (see Sec. |E 

The wave function of the original problem (|I|) can be obtained as a linear combination 
of ipn = \l/(— e~*^"g") where the momentum p_ belongs to a set of rational points of the 
Brillouin zone . For P = odd the rational points correspond to the center of a band. There 
are two distinct sets — one for t = +1, another for r = — 1. For P = even and r = 1 a 
rational point also corresponds to a center of the band. The rational points at r = — 1 and 
P = even, give bottom edges of odd bands and top edges of even bands for /i(— 1)-^/^ = — 1, 
counted from the bottom of the spectrum and top edges of odd bands and bottom edges 
of even bands for /i(— 1)^^^ = -|-1. This allows one to find bandwidths and gaps. Let us 
take P = even and numerate the energy levels corresponding to rational points with r = — 1 
according to their order starting from the bottom of the spectrum En, n = 1, . . . ,Q. Then 
the width of the n-th band (Wn) and the n-th gap (Dn) are: 



^We also restrict ourselves to the isotropic (critical) case A = 1. 
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Wr, = + ^Q+l_„), (13) 

^ \ Eq^n+i - Eq^n if P/2 + n - even ^ > 

For P-odd and Q-odd edges of bands and their widths may be obtained by the flux reflection 

q -q-\ 

The integrability manifests itself in the fact that the solutions to the difference equation 
(0) are polynomials 

i=0 

with roots obeying Bethe equations: 

Zi - qzj {q2Zi + iTK) {q2Zi- ik) 

These equations have Q independent solutions. Each of them consists of Q — 1 complex 
numbers {2:1, . . . , zq^i} and corresponds to an eigenstate of the Harper equation at a rational 
point of the Brillouin zone with the energy: 



E = ifiq^{q-q 



I- 



[l-q-Q){l-T) 



,=1 gi/2 _ q-1/2 



(16) 



The Bethe equations ( [T5| ) are similar to the equations for the XXZ chain with maximal 
allowed spin S = and commensurate anisotropy cos(^?7). The major difference is that 
it corresponds to a spin-problem on one or two sites [0,^]. 

Below we give evidence that in an incommensurate limit all Bethe roots are grouped to 
strings. Each string of the "spin" / is a complex of 2/ + 1 roots 

{zl!^ = vixiqr, m = -/,-/ + !....,/} (17) 

which have a common modulus X/ > (the center of the string), a parity Vi = ±1 and differ 
by multiples of qi. We show that qi is the 21 + 1-th root of unity, "closest" to q (Sec. |V|) 
and xi = 1 + 0{Q-^). 

Each state is characterized by a set of spins {lk,lk~i, ■ ■ ■} and parities {^^,^^-1, . . .} of 
the string content, such that the total number of roots I]f=i(2/i + I) = Q — I 

Ik 

^{ik,vk,...}^^^ Y[ (z-xiViq';^)^^^'-"'"'~''-^z). (18) 

m=—li. 

Moreover, in the next section we argue that the wave function with a "subtracted string" 
■qj{ik~i,vk-i,---} ^^'j is approximately the wave function of the AH problem of the parent gener- 
ation, i.e. the problem with flux r]' = r/(J~^) (^ and the length of the "subtracted" string is 
the Hall conductance Ik = \a{k)\. The formula ( p!7D is asymptotically exact with accuracy 
increasing with the length of the string as 1/(2/^ + 1)^. 
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Comparison with the XXZ chain may be useful. In that case strings become exact in 
the thermodynamical hmit, i.e. in the hmit where the size of the chain and number of spins 
up go to infinity faster than anisotropy angle P/Q approaches an irrational number. In 
that case spins (i.e. lengths) of the strings and their parities are drawn from a subset of 
integers, called Takahashi-Suzuki numbers, while centers are continuously distributed on 
the real axis. Our case is different — the number of sites on the equivalent chain is 2. Also 
one does not expect to describe strange sets by continuous distributions. Nevertheless, 
the string hypothesis remains robust. Naturally the incommensurability plays the role of 
"thermodynamical limit" . In the Sec. |V| we show that allowed spins of strings remain to 
be the Takahashi-Suzuki numbers, however the number of strings with a given spin is finite 
and their centers approach 1 with increasing length of the string. 

Strings play a bigger role than just being solutions of the Bethe equations. The spins 
of strings of which the state is composed give the topological characteristics of the state, 
related to the holonomy of the wave function in the Brillouin zone (the Hall conductance). 
They eventually describe the topology of the spectrum. 



IV. STRINGS AND HIERARCHY TREE 

In this section we present heuristic arguments for the string hypothesis and support them 
with numerics in Sec. |VI[ The arguments are based on the analysis of singularities of the 
Bethe equations (plSf). It turns out that the string solutions follow the hierarchical tree of 
Sec. f5 

Let us assume that a solution of the Bethe equations (|^) contains a long string of spin 
/ ^ 1 ( P^ ) . Let us substitute this string into (|TBp and drop all terms, which are not singular 



as qi q. Equations for the first and the last root of the string ViXiqi ' and ViXiq\ give 

KXiViqi\q - qi){. . .) ~ {xmqY^ - iKq^T){xiViqf^ + inq^) (19) 
nxiviqq\~^^{qi - g)(. . .) ~ {q^XiViq\ + inT){q^ XiViq\ - in). (20) 

where by (. . .) we denote the terms of the order one at qi — > q. Compatibility of these 
equations requires 

\qfq±l\ = 0{\q-qi\). (21) 

If qi = (^^P"IQ" is sufficiently close to g = e^'^^ i.e. \r] - P" /Q"\ ~ Eq. (|l]) may be 

satisfied only if 

qf+^ = ^1 (22) 

where =F corresponds to the sign in (|2TD. This gives Q" = 2Z + 1. We will refer to the strings 
with qf^^ = +1 as to "closed" strings as opposed to "open" ones with qf^^ = — 1. 

For open strings the equations (P^^PDj ) give an approximate value of the center ~ 1 
and two possible values of parity: 

Vi = -iq\^"\={-li'^'^K, (23) 
= Mi^^''^TK = -(-l)[''Vfi;. (24) 
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The symmetries of the Bethe equations require a center to be a real positive number. In 
case 7^ 1, there is always another string with an inverted center xf'^. 



If qf'^^ = 1 the leading singular terms ( ]T9| , |20| ) of the Bethe equations vanish. The center 



and the parity of the closed string are determined by the next leading singularities. We do 
not study these here. Empirically we still find that centers of closed strings are close to 1 
as for open strings. 

Now consider the Bethe equation for a root which does not belong to the string. Then, 
assuming qi to be very close to q after a set of cancelations, we obtain again the Bethe 
equations for Q — 21 — 2 remaining roots. If qf^^ = 1 the equations retain their form. In 
case qf"^^ = — 1 the Bethe equations for the point {t,k,ix) are transformed into ones for 
the point (— r, —k, fi) for a string with the parity (p3|) and to the point (— r, k, —fi) for the 
parity 

In other words, while subtracting the string from a wave function of a generation where 
rational points correspond to edges of bands (P is even) we obtain the Bethe equation for the 
generation where rational points are the centers of bands {P' is odd). However, subtracting 
a string from the generation where rational points are centers (P is odd) one may get P' 
both even and odd, depending on the sign in eq. (^2]) . 

Now comes the crucial step: we have to determine what values 2/ + 1, the length of 
the string and its parity can take. The new approximate Bethe equations on the remaining 
Q — 21 — 2 roots would be consistent if one replaces q by the closest q' such that (g')<3~2'-i = 
±1. This requirement and also the requirement that qi is also the closest to q, leaves us with 
two possibilities 



qi = e 1 qi = G '^-'^r 



p-Pr--^ or , p 



(25) 



q = e "^-"^r-i q = e 

This means that the states of generation P/ Q have the largest string of the length 

2/ + l = {g-g,_i or Qr-i} (26) 
The parent generation of a state may have a flux: 

respectively. Thus we see that the lengths of two possible highest strings are the same 
as two possible values of Hall conductance of bands (^. We conjecture (and checked this 
numerically) that the length of the highest string is equal to the Hall conductance of the 
band. 

Now let us discuss the parities of open strings. If r = — 1 the both formulae (p3| , |2^ ) give 
the same result v = {—ly^^'^K. However for r = +1 (center of bands) there is a choice between 
two possible values of parity. To make this choice we must further specify the hierarchical 
tree and the spectral flow. Let us consider for example the band k which belongs to the 
generation t] = P/Q with P-odd and r = +1. Then the corresponding parent band k' 
belongs to the generation r]' = P'/Q' with P'-even (we consider the case when the largest 
string is open) and r' = — r = — 1. There are two different energy levels corresponding to 
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the rational points in Brillouin zone which belong to the parent band k'. These are top and 
bottom of this band. We will call the state corresponding to the top (or bottom) of the band 
k' as parent state of the state in the middle of the k-th band if the Hall conductance in the 
gap adjacent to the top (bottom) of the band k' is the same as the one in the gap adjacent 
to the top (bottom) of the band k of the generation rj. If one goes from the generation rj' to 
generation rj in AH problem only new gaps open in the spectrum while the old ones present 
in parent generation stay and have the same Hall conductance in the daughter generation. 
Therefore among two states at the top and the bottom of the parent band we call the parent 
state the one which is not splitted of from the band by new gap while evolving from the 
generation rj' to the generation rj. The both states can be parent states if the band k' is not 
split by a new gap. 

The edge of the band k' corresponding to the rational point characterized by k' and n' is 



a top of the band if 



)P'/2+k'> 



jj,' = —1 and it the bottom if (— l)'^'/^"'"'''/i' = +1. Using the 



rules of transformation of k and /i while subtracting open string from the state in the middle 
of the band k one can easily obtain that the right choice of parity is vi = yUK(— 1)'^'/^+'^'+''''] 
if the parent state is the bottom of the band k' and vi = — /!/?(— 1)"^ Z^^*^ if it is the 
top. Finally, the k-th gap at the generation rj is an old one (present at the generation r]') if 
corresponding Hall conductance |crfc| < Q'/2. and we obtain: 



+ 1 
-1 



if 
if 



-i\<Q'/2 
<Q'/2 



(28) 



If both of formulae in (|28D work one can choose parity to be f = ±1. We notice than in this 
case the string content of the state in the middle of the band k must include the string of 
the largest length at least twice with both possible parities. 

Summing up we obtain the iterative procedure of constructing the string decomposition: 

• Starting from a band of the generation i]^^'^ = Pr/Qr and a point of the Brillouin zone 
T,K,,fi, we determine the string with highest spin (^). Its parity is given by (p3|). 
Subtracting this string we obtain the Bethe equation for a generation r]^^~^^ (p7D for a 



point of the Brillouin zone r' = qf^~^^T, k' = ±k and so on. Continuing this procedure 
to the end of the tree we obtain an unambiguous string decomposition of each level. 



p 

^,T, K,, fJ, 



the set of 



To illustrate the iterative hierarchical procedure, let us denote by 
roots {Q — 1 complex numbers) for flux P/ Q and rational points of a particular band labeled 
by r,K and fi. Then the iterative procedure can be symbolically written as 



Q 



{^^■(21+1)} 



P' 



-,qf-^'T,±K, ±/i 



Q 



(29) 



Here the symbol {vi ■ (21 +1)} denotes the string of the length 2/ + 1 with the parity vi, 
Q' = Q - 21 + 1, P' = P - P" and P" is the nearest integer to {21 + 1)1]. If the string is 
open, i.e. qf"^^ = — 1 the numerators P and P' have opposite parity (if P is odd, P' is even 

and vice versa), while P and P' have the same parity if the string is closed. The symbolic 

' p 



expression (p9D means that Bethe roots corresponding to the state 



are given by 
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the roots of the state 



plus the roots of the string of the length 21 + 1 
with parity vi. 

The iterative procedure of the string's decomposition is identical to the hierarchical tree 
based on the spectral flow (Sec. H). Once the length of the highest string in a given band is 
the Hall conductance of the band, the string length is the rate of the spectral flow from the 
parent bands. Subtracting this string we arrive at the parent band of the generation given 
by the second term of the sequence of intermediate fractions (^). Performing this procedure 
recursively, we obtain the entire set of possible lengths are the celebrated Takahashi-Suzuki 
numbers (|l3)) known from the Uq{Sl2) integrable models. 

Since lengths of strings are Hall conductances, the string content of a state can be read 
off from the hierarchical tree. The lengths of these strings are the Hall conductances of the 
bands on a branch (connecting these bands) that leads to the state. The string decomposition 
of the 15 bands of the example of the Sec. P/Q = 4/15 ([TT|) is, as shown in Fig. |l] (only 
the lower half of the tree is shown): 

-4), Jg = (±2,1,-7,-4), 
Jio = (±2,1,11), 
Jn = (±2,1,3,-4,-4), 
Ji2 = (1,1,4,-4,-4), 
Ji3 = (1,1,1,-7,-4), 
Ji4 = (1,1,1,11), 
Ji5 = (1,1,1,3,-4,-4), 



■h = 


(-1, 


-1,-1,-3,4, 


J2 = 


(-1, 


-1,-1,11), 


h = 


(-1, 


-1,-1,7, -4) 


J A = 


(-1, 


-1,-4,4,-4) 


J5 = 


(±2, 


-1-3,4, -4) 


Je = 


(±2, 


-1,11), 


J7 = 


(±2, 


-1,7, -4), 


Js = 


(±2, 


-4,4,-4), 



(30) 



where in brackets we showed the lengths of strings and their parities as signs in front of 
lengths. 

Fig. 1^ shows the root patterns of the 7-th branch (counted from the bottom of the 
spectrum) Jy of the tree from Fig. |l]. We added also the root pattern of the 15-th level 
for the flux 9/34 which belongs to the same sequence. One can notice that the subsequent 
patterns of roots differ only by the string determined by hierarchy procedure. Although the 
topology of the hierarchical tree seems to be exact, its quantitative description in terms of 
strings is approximate. The true values of roots deviate from their bare value xi = 1 (see 
Fig. H). Eq.(|l^) is approximate and valid only for long strings. 



V. GAP DISTRIBUTION AND BANDWIDTHS 

A. Gaps 

A direct application of the string hypothesis is the gap distribution p(-D), i.e. the number 
of gaps with magnitude between D and D + dD. This problem is essentially reduced to the 
scaling of the smallest gaps. Indeed, while moving along the tree new gaps appear, while 
older gaps become wider. There is no overlap between bands except at -E = for even 
denominators of the flux (they touch) — all gaps are open pT| . 

Thus while moving from a generation P' / Q' to generation Pj Q, there will he Q — Q' new 
gaps. For typically Diophantine flux denominator grows fast with the generation number 
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and Q — Q' ^ Q. This means that the integral J^™'"' p{D)dD = Q is saturated by the lower 
limit. If p ~ D~^, then D ~ Q. 

To find the size of the minimal gap, let us consider P-even and r = — 1. Then the Bethe 
equations give edges of gaps. Let us consider an arbitrary gap in this generation and trace 
the genealogy of lower and upper edges of the gap along two paths J_ and J+ of the tree, 
until we find the nearest common ancestor at the generation with denominator Q' . The 
minimal gap would correspond to the shortest paths J±, i.e. when Q' is the denominator of 
the generation parent to P/Q. The string decompositions of the two edges have a common 
part consisting of strings of length smaller than Q' (they correspond to the common path 
from the origin of the tree to the parent generation Q') and different strings with length 
bigger than Q' belonging to the paths J±. Then according to (|1^), the width of the gap is 
the difference between energies of strings with length greater than Q' along the path J_ and 
the energies of strings along J+. The contribution to the energy ei from a string with a spin 
/ which connects the center of the band (r = +1) with the edge of the band (r = — 1) is 

where we used the fact that for a typical Diophantine flux 77 the length of the largest string 
21 + 1 scales as Q along the sequence of rational approximants. This gives us an estimate 
for the gap Ei ~ The difference of energies of strings is of the same order. This 

gives -Dmin ~ l/*?^ aiid the gap distribution function p(l^) ~ V^^/'^. This result has been 
obtained numerically in Ref. 



B. Bands. Stating a Problem 

The bandwidths are a more subtle matter. Let us take k-ih band belonging to some 
generation with P-even at r = — 1. Then is the edge of this band. We can take one step 
back along the hierarchical tree. We flnd that the parent band k' is of generation rj' = P'/ Q' 
with P'-odd and r = +1. The rational points at this generation give centers of bands. If we 
had started from the band Q + 1 — k (symmetric to the k-th band) we would flnd the band 
Q' + 1 — k' as its parent band. The roots of symmetric bands of generation r]' are related 
by symmetry Sec. IXp. Thus the string content of two paths leading to the 



bottom and the top edges of the band differ only by the highest string of the same length 
/ and the same parity vi. If we assume that adding a string does not change the centers 
of already existing strings (the noninteracting strings approximation), then the bandwidth 



is given just by twice the energy of the highest string (|3ll) . This gives an estimate of the 
bandwidth ~ Corrections to this crude estimate come from many sources (roots get 
some small corrections when the new string is added, every root has "flnite size" corrections) 
and these are of the order as well. 



It has been shown by Thouless |33] that the sum of all bandwidths scales as (more- 
over \imQ^ooQJ2k^k = const = 9.32...). If all bands scale similarly (fractal), then the 
scaling dimension would flt the scaling of the total bandwidth. We call this the mean di- 
mension. In reality, bandwidths scale differently (multifractal). Their anomalous dimensions 
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ej depend on the path J of the hierarchical tree that the band belongs to: W{J) ~ Q 
where Q is the denominator of the flux running along the path. According to ref. |^ for 
?7 = exponents ej vary between 0.171 and —0.374 , where the positive anomalous di- 
mension corresponds to a path in the center of the tree, while the negative one belongs to 
the path along the lowest bands. 

A distribution of anomalous dimensions of the bandwidths, can be described by the 

function: = — limQ_,oo °^ i^gQ — ~ Here is a limited list of what we do know about this 

function. The function is convex and has linear asymptotes = for /5 —oo 

and F{(3) = for (3 oo. A deviation of F{(3) from a straight line (fractal) indicates 
multifractality. The point f3 = Z)h where F{D^) = is known as the "Hausdorff dimension" 
and the slope of F at this point, gives the most probable scaling dimension eo = -F'(Dh)- 
The scaling of the total bandwidth |To|j33[] gives F{—1) = 1. Obviously, F{0) = —1. We 
have plotted F{i3) for the case where the flux is equal to the golden mean in Fig. ^ 

If F{/3) were a straight line, Dh would be | and eg = 2. It would mean that each scaling 
dimension e(J) is equal to the mean field value 2. From Fig. |^ we can see that F[f3) is 
indeed close to a straight line. The Hausdorff dimension = 0.49792 ± 4 x 10~^ and the 
most probable scaling dimension eo = 1.9996 ± 1 x 10~^ are numerically close to | and 2 
respectively. Anomalous dimensions are small but code the most interesting features of the 
problem. 

The topological classification of solutions to the Bethe equations by means of strings 
alone gives the mean-field dimension of bandwidths and absolute values of bandwidths with 
high accuracy of the order of 1/Q^. Analytical determination of anomalous dimensions goes 
beyond just the topological classification and is perhaps the most challenging problem in 
the field. We do not attempt to solve it in this paper. 



VI. NUMERICAL EVIDENCE FOR THE STRING-HIERARCHY HYPOTHESIS 

In this Section we present a fragment of an extensive numerical analysis we have per- 
formed in order to confirm the string-hierarchy hypothesis of the Sec. |V[ We choose the 
most "irrational" flux (golden mean) rj = = [1; 1, 1, 1, 1 . . .]. This is the only case where 
intermediate fractions are the same as convergents 

Pi _ Fi^i _ 1 1 2 3 5 8 , . 

g; ~ "fT ~ T'2'3'5'8'T3''"' ^ ^ 

where Fi are Fibonacci numbers as well. They give the fluxes of generations. The set of 
Hall conductances = Takahashi- Suzuki numbers = allowed lengths of strings and are all 
Fibonacci numbers: 

2/fc + 1 = Qk-i = Fk-i = 1, 1, 2, 3, 5, 8, ... . (33) 

As an illustration let us consider the bottom edge of the lowest energy band of the 
generation ||. According to the string hierarchy hypothesis the length of the largest string 
is 21 and then iteratively one obtains that the string content consists of strings of lengths 
-^3fc+i = 1,5,21 for /c = 0,1,2 and can be written as (1,1,5,5,21,21). All string parities 
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in this example are —1. The actual numerically calculated roots are shown in Fig. ^. One 
notices immediately that the roots organize themselves into distinct groups whose members 
have roughly the same radius: two groups with 1 root(+), two with 5 roots(x), and two 
groups with 21 roots(*). 

The spins of strings and their parities are topological characteristics of the state. In 
contrast to these properties, the centers of strings are not topological and become close to 
1 only for large spin strings. For our example they are x = 2.10, (2.10)"^, 1.23, (1.23)~^, 
0.81, (0.81)~^, 1.05, (1.05)~^. We have found numerically that the centers of large strings 
are xp^ ~ 1 + tt^ — with f = 1.14. 

Another example. Fig. ^, illustrates the hierarchy. It shows how 12 roots of the band 
edges of the generations ^ can be approximately obtained from 4 roots of the centers of 
the bands for the parent generation | and 7 roots of generation by the addition of the 
strings of the lengths 8 and 5 (□) respectively. 

The string decomposition gives an asymptotic form of the wave function (|1^). In partic- 
ular the wave function of the bottom edge of the lowest band of the generation Fs^.i/Fs^ 
is 

^{F3.-.,...}(^) ^ n - ^F.,.,,qs.^''^){z - x-^l^^^q^,^^'^)^^^^''---^ (34) 

i=k m=l 



where qj = e . This formula is asymptotically exact at A; ^ oo, but also works very well 
even for k = 0, ip^^^ = 1 and ipi^3k-2,---} jg ^}^g ^ave function of the generation -Fsfc-iZ-Fsfc. 

First of all let us check that the formula (plj) is qualitatively correct and gives the 
same root pattern as the actual roots of the Bethe Ansatz equations. In Fig. ^ we have 
presented, in addition to the numerically calculated roots corresponding to the ground state 
of equation (|T2p for the flux r] = 34/55, also the roots given by (|5^). It can be seen that the 
roots produced by ( ^41) are in quite good agreement with the numerical results. Actually 
one can check that the numerical accuracy with which formula ( ^4|) gives roots is ~ 1//^ for 
each root of the string of the length /. Now let us calculate the overlap between this wave 
function (|3^ ) and the direct numerical solution of the eq.(H)- 

The result is tabulated in the following table, where the overlap between the wave func- 
tion obtained from the exact (numerically obtained) roots V'Exact and ansatz wave function 
of different order k: V'Ansatz shown. The ansatz wave function of the order k is given by 
eq. (p^ . It is "made" out of real (numerical) and ansatz roots. The order k of the ansatz 
function tells how many of the longest strings are replaced by ansatz ones, k = corresponds 
to pure ansatz wave function. The bigger k the more short strings are replaced by ansatz 
strings. The positions of centers of the ansatz strings are taken to be xp^i+i = 1 + ^/F^i+i 
with the only fitting parameter ^ which is chosen to be ^ = 1.14. 



k 





1 


2 


3 


4 


5 


(^Ansatz I^Exact) 


0.99475 


0.99464 


0.99657 


0.99669 


0.999435 


1.000000 



The last column of this table reflects trivially the normalization of wave function. One 
can notice that for A; = 4 when two strings of the length 377 are replaced by the ansatz roots 
the overlap between ansatz wave function and exact one is 0.9994 which is very close to 1. 
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Actually for this case even the purely ansatz wave function (with all roots produced from 
(H)) has a good overlap: 0.995 with the exact one. 

A more subtle check of the high degree of similarity between the exact and ansatz wave 
function is the distribution of scaling exponents of moments of the wave function, which is 
often used to characterize the multifractal properties of the set generated by of values of ipn 



35| . This distribution is given by the Legendre transformation 



of the function 



F(f3) = hm ^°g^"=il^"l^ (36) 

Q^oo logQ ^ ^ 

An extended support of /(a) gives some measure of multifractality. 

In Fig. ^ we have plotted /(a) for ^/^ Ansatz from eq. ( plj) and ^/'Exact from the exact roots, 
calculated numerically for a particular value of the flux, rj = 610/987. By representing the 
data in this way we have neglected taking the limit in eq. (p6|). For this value of the flux 
this is justifiable by the fact that when increasing the denominator of the flux beyond 987, 
F{f3) does not vary visibly. The figure shows that the ansatz approaches the exact f{a) 
remarkably well, with just one free parameter ^ ^ 1.14. 

Concerning the support of /(a), ttmin is approximately 0.32 for the ansatz and 0.31 for 
the exact wave function, whereas Omax is approx. 3.5 for both ansatz and exact wave function 
within margin of error. 



VII. CONCLUSIONS 

The Azbel-Hofstadter problem as a typical quasiperiodic system generates a complex 
spectrum. Similar complex spectra were found in a number of dynamical systems. Since 
the empirical observations of Hofstadter 0, evidence has been mounting that these spec- 
tra are regular and universal rather than erratic or "chaotic". They have a deterministic 
hierarchical structure. In this paper we have shown that the topology of the strange sets, 
generated in the problem is determined by the Chern numbers of the spectral curve i.e. by 
the Hall conductances. Even more so — at every finite step of the hierarchy the spectrum is 
integrable. We have found an anticipated match between Hall conductances and dimensions 
of representations of the quantum group Uq{sl2)- The latter are Takahashi-Suzuki numbers 
or the lengths of the strings of the Bethe Ansatz solution. This correspondence suggests a 
natural hierarchical tree, which, we believe, is relevant for general quasiperiodic systems. 

In the Bethe Ansatz language, each state is characterized by a particular string con- 
tent. Proceeding along the tree toward the incommensurate limit corresponds to addition of 
strings. This picture is somewhat reminiscent of the discrete renormalization group approach 

We were concentrated on the topological aspect of the string hypothesis. It alone gives 
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• The explicit asymptotically exact form of some wave functions for irrational flux. 
Multi-fractal properties of these functions, although not exact, are in good agreement 
with numerical results. 

• The scaling exponent 3/2 of the gap distribution. 

• The most probable scaling dimension 1/2 of the bandwidths. 

However, the major ends of the strings are loose. The string hypothesis solves the Bethe 
Ansatz equations with an accuracy 0{Q~'^), i.e. is asymptotically exact in the incommen- 
surate limit Q ^ oo. However, the most interesting quantitative characteristics of the 
spectrum are actually in the finite size corrections of the order of to the bare value 

of strings. Among them are the anomalous dimensions of the spectrum. We believe that it is 
possible to find them analytically via a more detailed study of the Bethe Ansatz equations, 
beyond just the analysis of singularities. This is a technically involved but a fascinating 
problem. The ultimate solution of the problem, however, would be through the applica- 
tion of conformal field theory approach, which has been proven to be effective for finding 
the finite size corrections of integrable systems, without the actual solution of the Bethe 
Ansatz. 
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IX. APPENDIX: BETHE- ANSATZ EQUATIONS FROM HOFSTADTER MODEL. 



In this section we (re) derive the Bethe Ansatz equations for the Azbel-Hofstadter problem 
with an anisotropy parameter A |28|,^ . 



The Hamiltonian of a particle on a lattice in a magnetic field can be written as 

H = T, + T-' + XiTy + T;'), (37) 

where the translations operators on a square lattice in a magnetic field 2711] per plaquette 
form a Weyl pair 

T,Ty = q%T,. (38) 

The representation T^fn = fn+i, Tyfn = q^'^e^^y fn with q = e™'^ gives the Harper equation 
(0). If the flux through a plaquette of the square lattice is a rational fraction of the flux 
quantum r} = P/Q then the operator commutes with the Hamiltonian ( pT]) and there 
are finite dimensional representations characterized by two momenta k^, ky so that 
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J^Q — ^iQkx rpQ _ ^iQky ^gg^ 

The non-equivalent representations are given by restricting < k^. < 2tt/Q and Q < ky < 2% 
and Schrodinger equation becomes: 

Eijn = e'^^i)n+i + e-'^^i)n-i + 2A cos(27rr/n + ky) (40) 



with periodic boundary conditions V'n+Q = V'n- The Chambers relation for the problem 
(^) says that the energy spectrum is obtained as the solutions of an algebraic equation 
Pol{E) = A{kx, ky) where Pol{E) = + . . . is a polynomial of Q-th degree which coeffi- 
cients do not depend on fc^ and ky. The only dependence of energy on k^ and ky comes from 
free term of algebraic equation via the combination: 

A{k^, ky) = 2 cos Qk^ + cos Qky. (41) 

Therefore the edges of the energy bands are given by extrema of A which assumes a mini- 
mum/maximum given by A = ±2(1 + A*^). 

Below we use another representation (the modified chiral gauge): 

U -\- a 1 1/ + a , ^ 

T. = UV-^, Ty = VU-'j^, UV = qVU. (42) 

In this representation for the choice 

a = —iKA 2q 2^ 
b = —iTK\2q~2^ 



1 1 , 



where t, k = ±1 the Hamiltonian (|37| ) becomes 

H = qVU~^{U - iKX~^q"^){U + irnX^q"^' 

+Xq-^V-^U-\U + iK\-^q^){U - iTKX^qi). (43) 

Let us choose the following representation of operators U and V: 

Then the Schrodinger equation corresponding to (^) becomes 

- nq^'e-'P-X-^E^n = iqe'P+iq'^e-'P- - ir/tA-^g-^)(g"e-*P- + mA^g-^)z^„+i 

- iq-^e-'P+iq'^e-'P- + iTKX-^q^{q''e-'^~ - mA5g5)^„_^ (44) 



The equation (|4^) has periodic coefficients and one can require that i'n+Q = ipn i^n+2Q = 4'n) 
for P-even (P-odd). Then the solutions will be labeled by < p+ < 27r/(5 for P-even and 
by < p+ < tt/Q for P-odd. 

Let us now extend the periodic function ipn to the entire complex plane by writing 
tjjn = ^{z)\^^_^-ip_g„. We obtain: 
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iq ^P+{z — iTK,X ■^q'^){z + iK,X^^q^^)'^{q ^z). 



(45) 



If e^*P+ = 1, the eq. ( ^5]) has Q polynomial solution of the order Q — 1. Below we denote 
/i = e*^+ = ±1. The change /i — /i changes the sign of energy corresponding to given 
polynomial solution of (^). Let us parameterize the polynomial solutions by its roots zf. 



(46) 



i=l 



Evaluating eq. 



at one of the roots Zi, we obtain the Bethe Ansatz equations: 



Q qZi-Zj _ [zi-iTKqlX 2) (^z, + iKqlX2^ 

^ 11 / 1 '. 1 \ / 1 '. TTV 

j=i Zi ~ Q^j iq2Zi + tTK,X 2 j lq2 Zi — tK,X2 j 



(47) 



These equations have Q independent solutions. Each of them consists of Q — 1 complex 
numbers z^, k = 1, . . . ,Q — 1. Comparing the coefficients of z'^ and z on the left and right 
sides of equation ( ^5]) one finds the energy: 



E = ifiq^{q-q-^)X^/^ 



Q-i 

i=l 



qi/2 _ q-i/ 



(48) 



or 



E = in(q — q )X 



Q-1 

i=l 



gl/2 _ g-l/: 



(49) 



Combining these two formula we have a more convenient form of the energy 

0-1 



E = nx-q'i{q-q-')X'l^ 



i=l 



gl/2 _ q-1/2 



yl/2 _ ^^-1/2 



(50) 



The Bethe equations (|47D possess some discrete symmetries: if the set {zi} is a solution 
with anisotropy A and energy E then {z*} is also a solution with the same energy and 
{— rZj"^} is a solution with the energy q'^E. In its turn the set {z~^} is a solution of the 
problem with A^^, its energy is —rq'^E. Finally, there is a trivial symmetry of equations 
Zi, K —>■ —Zi, —K, which does not change the energy. 

After some algebra one finds a relation between the wave function of the Harper equation 
(PH) and polynomial (^). 

(i) In the case of g*^ = — 1 (P-odd) and in the middle of bands the relation is 



n —n2+2n 



-lYq 



2Q-lm-l 

E n 

m=0 j=0 



^ikyq2n+ 2 



1 _ k^ + ky_ _ 1 

1 A2 — TKe ' n J 



<1 



1 ■kx+ky_ 1 

A 2 — Ke 2 qJ^2 



^(-ie'^g™), (51) 



where k^, ky satisfy the condition A{kx, ky) = 0. 
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(ii) In the case of (p — +1 (P-even, Q-odd) and A = 2//(— 1) 2 (A*^ — r) the relation is 

0-\m-\ / \i • _,-^£+^ -1-1 \ 

V'n = g " EH e''^g2"+2 ^ |^(-e' ^ g"^)- (52) 

m=o j=o \ Aa+iree' 2 y 

If T = +1 it corresponds to the middle of bands and the rational points kx, ky given by 
the equation "'t|^^"f,^ = fi;(-l)^A2 . If r = -1, A = 2//(-l)l (A« + 1) the relation 

sin -2"(fes/+''i'2Q) 

p / ^ \ 

corresponds to the points at the edges of bands k,j. = ^ ^~^~2^^^ and ky = ^ ( _^ 2/ j 

for / = 0,1,...,Q — 1. If counted from the bottom of the spectrum, these edges are ordered 
as bottom-top-bottom. . . if 1)2" = —1 and top-bottom-top. . . if 1) 2" = 
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FIGURES 



Ti = 0/1 Ti = l/2 Ti = l/3 Ti = l/4 Ti = 2/7 Ti = 3/ll Ti=4/15 

T = -l x = +l t: = +1 t = +1 x = -l x = +l T = -l 

K = -l K = +l K = +l K = +l K = -l K = +l K = -l 




FIG. 1. The lower half of the hierarchical tree for the flux rj = 4/15 is shown. Nodes which lie 
on a vertical slice of the tree are bands of a generation ordered with respect to their energies. The 
numbers in the figure are Hall conductances of bands and gaps. The Hall conductances in gaps are 
shown as encircled numbers. "Top" and "bottom" are the edges of a band. For all states shown 
= +1. Root patterns of the highlighted branch J7 of the tree are shown in Fig. ^. 
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ri = 9 / 34 
x = +l 

K = +l 

Band= 15 



Ti = 3/ 11 
x = +l 

K = +l 

Band = 5 



ri = l/3 
x = +l 

K = +l 

Band = 2 




FIG. 2. The root patterns along one branch ( J7) of the hierarchical tree of Fig. ||. Subsequent 
root patterns along the tree differ by a string whose length is equal to the Hall conductances of the 
"younger" -one of the two bands. The parities of these strings are in agreement with (23) and are 
shown as signs in front of string lengths. The band numbers are given as counted from the bottom 
of the spectrum. For all states shown fi = +1. 
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FIG. 3. Shown in the figure are the roots that parameterize the ground state of the Hofstadter 
problem at flux 34/55. The roots are organized in strings: two of the length 1(+), two of the 
length 5(x) and two of the length 21(*). These roots can be approximated by ansatz ( |3^ and are 
plotted in the Figure as boxes(Q). Strings or more generally roots tend to repel each other. This 
"splitting" is captured empirically by the expression for Xn in (p^). 
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FIG. 4. This figure shows the evohition of root's pattern along the hierarchical tree: root 
patterns of the generation with the flux ^ are constructed out of previous generations with the 
flux I and |) by adding strings with the lengths 5 or 8 respectively with parities it denoted as sign 
in front of string lengths. The symbol ^, —0+ denotes a state at the bottom (— ), middle (0) or 
top (+) of the i-th band (counted from the bottom of the spectrum) of the generation with a flux 
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FIG. 5. The function f{a) for the wave function of the bottom of the lowest band. The sohd 
hne is an exact result (found numerically), the dashed line is the f{a) of the approximate wave 
function (|34[). The support of f{a) is 0.3 amm < a < Omax ~ 3.5. 
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FIG. 6. The free energy F{/3) is plotted as a function of /3 for the flux equal to the golden mean 
rj = The deviation from a straight line indicates multifractality. The leftmost inset zooms 

into the region f3 G [0, 1]. It includes the points F{P = 0) = — 1 and F{p = 1) = 1 known exactly. 
The rightmost inset zooms into the region where F{(3 = Dn) = 0. Dh is the Hausdorff dimension. 
The slope of F at this point gives the most probable scaling dimension ao = jt^^- Although both 
Hausdorff dimensions and the most probable scaling dimension are numerically close to 0.5 they 
differ from 1/2 within the numerical error-margin. The error margin (10~^— 10~^) is not visible in 
the figure. 
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